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Abstract. The article discusses the following frequently arising question on the 
spectral structure of periodic operators of mathematical physics (e.g., Schrodinger, 
Maxwell, waveguide operators, etc.). Is it true that one can obtain the correct spectrum 
by using the values of the quasimomentum running over the boundary of the (reduced) 
Brillouin zone only, rather than the whole zone? Or, do the edges of the spectrum occur 
necessarily at the set of "corner" high symmetry points? This is known to be true in 
ID, while no apparent reasons exist for this to be happening in higher dimensions. In 
many practical cases, though, this appears to be correct, which sometimes leads to the 
claims that this is always true. There seems to be no definite answer in the literature, 
and one encounters different opinions about this problem in the community. 

In this paper, starting with simple discrete graph operators, we construct a variety 
of convincing multiply-periodic examples showing that the spectral edges might occur 
deeply inside the Brillouin zone. On the other hand, it is also shown that in a "generic" 
case, the situation of spectral edges appearing at high symmetry points is stable under 
small perturbations. This explains to some degree why in many (maybe even most) 
practical cases the statement still holds. 
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1. Introduction 

The article discusses the following frequently arising question on the spectral structure 
of periodic operators of mathematical physics, which in particular is prominent due to 
the recent surge in studying photonic crystals [l9]-[2T], [36j-[38]. Let us have a periodic 
elliptic self-adjoint operator L{x,D) (e.g., Schrodinger, Maxwell), where we use the 
standard notation D = jV. The operator is considered in the whole space M", or 
in a periodic domain (on a periodic periodic manifold), e.g. in a periodic waveguide. 
The standard Floquet-Bloch theory (e.g., [H |35l |53]) shows that the spectrum of L 
in the infinite periodic medium can be obtained as follows: one fixes a value k of 
the quasimomentum in the first Brillouin zone B, finds the (discrete) spectrum of the 
corresponding Bloch Hamiltonian L{k) = L{x, D + k) acting on periodic functions, and 
then takes the union over all quasimomenta in the Brillouin zone. The question we 
address in this work is whether the correct spectrum can be obtained as the union over 
the boundary of the Brillouin zone 

This is well known to be true in ID (e.g., [131 ESj). In particular, the edges of 
the spectrum occur at the spectra of the periodic and anti-periodic problems on the 
single period. If this claim is correct in higher dimensions, the computational task 
is significantly simplified, due to reduced dimension. This is important, for instance, 
in optimization procedures, when one needs to run the spectral computation at each 
iteration [9l [10]. An experimental observation is that in most practical cases this 
is correct. One frequently encounters the belief that this is always true (albeit no 
justification is ever provided). On the other hand, unlike in ID, there is no analytic 
reason for this property to hold. Moreover, many researchers are aware that numerics 
sometimes produces counterexamples. Surprisingly, such examples are hard to come 
by and are usually not very convincing for an analyst (e.g., the error in computing 
the spectrum using only the boundary of the Brillouin zone is usually very small). 
The experience is that one needs to make the medium inside the fundamental domain 
(Wiegner-Seitz cell) truly asymmetric to achieve such examples. 

The first goal of this text is to provide simple definite examples to disprove the claim 
that the edges of the spectral bands can be found by using the boundary of the (reduced) 
Brillouin zone only. This is done by first analyzing some discrete graph systems. Section 
|2]describes such combinatorial graph counterexamples. Sect ion [3] deduces from this some 
quantum graph (see [39]) examples. Then in Section HI we bootstrap this to examples of 
waveguide systems or Laplace-Beltrami operators on thin tubular branching manifolds. 
Possibilities for obtaining counterexamples of the Schrodinger and Maxwell cases are 
discussed in Section |5l 

I If additional symmetries are present in the system, one considers the reduced (with respect to these 
symmetries) Brihouin zone. Another version of this question is whether the edges of the spectrum are 
attained on the set of high symmetry points of the BriUouin zone only. The importance of such points 
has been known since at least the paper When one needs to find the density of states, the full 
Brillouin zone is always required. 
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It is surprising, however, that the claim that we show to be incorrect in general, is 
still correct (or almost correct) so often. Thus, in many practical cases, computations 
along the boundary of the (reduced) Brillouin zone (and as a matter of fact, often at high 
symmetry corner points only) provide the correct spectrum. We suggest an explanation 
of this effect in the final Section [6l There we attempt to explain how it can happen that 
one often sees the spectral edges occurring at the high symmetry points only. It is shown 
that "generically" this occurrence is stable under small perturbations. In other words, 
there are open sets of "good" and "bad" periodic operators, the boundary between 
which consists of non-generic operators. This probably explains the frequent occurrence 
of the effect in practice. 

Finally, the last sections provide additional remarks and acknowledgments. 

2. Combinatorial graph examples 

We start by considering difference operators acting on a periodic graph. These will 
serve to illustrate the general ideas in a situation which is not difficult to analyze. 
Furthermore, building upon them, we will provide examples of more complex periodic 
spectral problems with the desired spectral feature. 

2.1. The main graph operators 

We consider the Z^-periodic planar graph F, with the fundamental domain W shown 
in Figure [H below. One imagines the graph F as obtained by tiling the plane with the 

#6 




is 

Figure 1. The graph F with fundamental region W. 



Z^-shifted copies of W. We will label the vertices in W and near W with the numbers 
shown in Figure [H 

Let ^^(F) be the Hilbert space of square-summable functions defined on the set 
of vertices of F. The discrete Laplacian on F can be defined in several (not always 
equivalent) ways (e.g., [6, 7J). We will use two of these. 

The first one. A, is defined for / G ^^(F) by 

(A/)(t;):=$^/(«), (2.1) 
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where the notation m ~ t> means that vertex u is adjacent (connected by an edge) to 
vertex v. For instance, 

(A/)(5) = /(l) + /(2) + /(4). 

The Laplacian defined in (12. ip differs from another discrete Laplacian often used in 
the hterature by the term d^fiv), where d^ is the degree (valence) of the vertex v: 

We will also employ another version, L, of the Laplacian, which is defined as 

(^/)(^):=^E^/(^)- (2.2) 



dn; \/ dni 

One could call it the Laplace-Beltrami operator. The need for this operator in our study 
will become clear when we move to the quantum graph case. 

One notices that both A and L are bounded operators in ^^(F). 

In the case of a regular graph (i.e. the one with constant degrees d^ = dof vertices), 
the spectra of A and L can be easily related. However, our graph F is not regular, and 
thus these spectra need to be studied independently. 

The following statement is well known (e.g., [6]) and immediate: 

Lemma 2.1. The operators A and L commute with any automorphism T G Aut(F) of 
the graph F. In particular, they commute with l?-shifts on F. 

For p = {pi,P2) G Z^, we denote by T{p) the shift operator by p on F. E.g., T(l, 0) 
shifts the vertex 4 to the vertex 7 and T(0, 1) shifts 3 to 6. 

Due to the periodicity of the operators, one can use the standard Floquet-Bloch 
theory [H [T31 ES] to study their spectra. In the particular case of graphs, this theory 
is also described, for instance, in [16], [3l]-[l0], ^9] . 

Let k = (fci, be a quasimomentum in the Brillouin zone B = [— tt, tt]^. Consider 
the space £| of all functions / satisfying the following Floquet (Bloch, cyclic) condition: 

f{T{p)v) = e'^-'f{v) (2.3) 

for all p G Z^. Here p-k = piki +P2^2- Such a function / is clearly uniquely determined 
by the vector (/i, /2, /s)* of its values at the five vertices in W , and thus £^ is five- 
dimensional and naturally isomorphic to £^{W). 

Definition 2.2. We will denote by H the boundary dB of the Brillouin zone B = 
[— 7r,7r]^ and by X the set of points k = (/ci,/c2) £ B such that ki and k2 are integer 
multiples ofn. 

We now define the Floquet Laplacian A(A;) : i'^{W) — * i'^{W) as the restriction to 
the space ij. of the operator A defined as in (12. ip . In terms of the basis of the delta 
functions at vertices of W, this operator has the following matrix: 



A(A:) := 
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In a similar way, one can define L{k) and observe that 
L{k) = S-^A{k)S-\ 



(2.5) 



where S is the matrix 



S:= 


V 





V3 














\ 





J 



[2.6) 



We now state the standard conclusion of the Floquet theory about the relation 
between the spectra of A and A(A;), or Land L(fc) (e.g., [l3l[T6l[3l[35lll0lll9l[53]). For 
each fixed k, the matrix A{k) (correspondingly, L{k)) is self-adjoint, and thus admits a 
spectrum of 5 eigenvalues {^j{k)}^^-^ (correspondingly, {/ij(fc)}^=i), which we number 
in non- decreasing order. It is well known that then each of the functions Xj{k) is 
continuous. The multiple- valued function k i— > {Xj{k)} is called the dispersion relation. 
Its graph is the dispersion curve, also called the Block variety. Each of the individual 
functions k Xj{k) is called the jth branch of the dispersion relation. 

Proposition 2.3. /Tgj. ISR The spectrum of A (correspondingly, L) is the 

union over k E B of the spectra of A{k) (correspondingly, L{k)): 

a{A)=[ja{A{k))=[j[j\,{k), 
a{L)=\Ja{m)=\J[jf^,{k). 

keB keB j=l 

The segments Ij = [J \j{k) (and their analogs for the operator L) are called bands of 

keB 

the spectrum of A (correspondingly, of L). 

Notice that our graph F does not have any point symmetries, and thus the reduced 
Brillouin zone is equal to the full one. So, the question we would like to address is 
whether one can replace the union over A; G -B in (12. 7p by the union along the boundary 
S = dB of the Brillouin zone B only. A more restricted question is whether the band 
edges are attained at points of X only. As we will show in the next sub-section, both 
of these properties do not hold, and calculations along the boundary lead to significant 
errors in spectra of A and L. 



2.2. Spectral edges - counterexamples 

In this sub-section we show that computations along S = dB (and thus over X 
as well) do not necessarily lead to the correct spectra of A and L, and the errors 
can be significant. We are interested in whether the segments Ij = IJfceB'^i(^) 
— UkedB^jW coincide. We will show that for our examples, even the unions 
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5 5 

o"(A) = IJ Ij and IJ Jj are different (the situation is analogous for L). This means 

that using only the boundary of the Brillouin zone, one does not recover the spectral 
edges (and thus the spectrum) correctly. 

2.2.1. Operator A One can easily find the spectrum of the Floquet Laplacian A(A;) 
(which is a simple 5x5 matrix), using for instance Matlab. Computing it for a grid 
in the Brillouin zone (we have used the uniform 64 x 64 grid in 5), one can obtain the 
whole spectrum of A. This leads to the following numerical values of the five bands: 
[-2.73,-1.90], [-1.63,-1.00], [-0.73,0.73], [0,1.46], and [2.00,3.23]. One notices that 
there are spectral gaps present between all consecutive bands, except the 3rd and 4th 
ones, which overlap. 

Since it is sufficient for our purpose to provide a single counterexample, we will 
focus on the second band [—1.63, —1.00] only. 

A gray scale plot of the second branch (corresponding to the second band of the 
spectrum) is given in Figure [2][§|. This numerical evidence shows that the band edges 



-1.001 



-1.630 



Figure 2. A gray scale image of tlie second brancli of tlie spectrum of A(fc). Extrema 
points are higliliglited. 



(i.e., the extrema of the branch function) occur at some values k not in S. This is 
confirmed by the graph of the branch presented in Figure [31 Let us now make this 
observation rigorous by finding the maximum and minimum values of A on S. The 
characteristic polynomial of A(/c) is 

ca(A; A;) = A^ - 8A^ - (2cosA;i + 4cosfc2 + 2)\^ 

+ (8 — 2 cos(A;i + A;2) — 4 cos ki — 2 cos k2)\ 

-2 cos(A;i + k2) -2 cos(A;i - A;2) + 4 cos k2. (2.8) 

Since the second band does not intersect any others, standard perturbation theory [32] 
implies that the corresponding eigenvalue branch A(A;) = \2{k) is analytic. It is not hard 

§ This and other plots are drawn over the square [0, 27r]^, rather than the Brillouin zone B — [— tt, tt]^. 
The origin (0, 0) is located in the upper left corner. 
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Figure 3. A graph of the second branch for the graph F. 



to check all possible extremal values of A on S. Indeed, assuming that ki equals ivr or 0, 
one can differentiate the secular equation ca(A; k) = with respect to ^2 and use that at 
an extremal point (unless it is one of the points of X), one has ^ = 0. One can do the 
same with the roles of ki and /c2 reversed. This calculation shows that extremal points 
can be located only where ki and /c2 are integer multiples of vr, i.e. at X. All points of 
X can be checked to yield that the minimum value on X is — \/2 ^ —1.414, attained at 
k = (tt, 0) and k = (tt, tt), and the maximum value is —a/4 — y/S ~ —1.082 at A; = (0, vr). 
These values can be compared with the numerically found extreme values of —1.630 at 
k ^ (1.865,0.785) and -1.001 at A; ^ (-1.080,5.203). This gives a difference of about 
8% at the upper edge, and 15% at the lower edge. These values and their symmetric 
counterparts are highlighted in Figure [21 

In fact, the location and value of the maximum of the second branch can be found 

— , — j is a 

maximum. 

Thus, we have an example of the situation when restricting the search of the edges 
of the spectrum to quasimomenta from H, leads to significant errors. 

One can draw other branches of the dispersion relation. They show that some of 
the branches do attain their extrema on X only, while some others do not. 
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2.2.2. Operator L The operator L can be analyzed similarly. We briefly summarize 
the findings. The characteristic polynomial for L{k) is 

1 



cl(A; k) 



X' - -X' - 
9 



- + icosA:, + lcosA:, 



+ 



13 



7 11 

cos ki cos ko cosf/ci + ko' 

162 162 54 ^ 54 ^ ' 



) A 



+ — COS /C2 
oi 



1 



cos(A;i — k2 



cos(fci + k2). 



(2.9) 



162 ' ' 162 

We observe that fl2.9l) is quite similar to fl2.8l) . modulo changes in some constants. The 
structure of the branches of solutions to cl{X; k) = are also qualitatively similar 
to the ones for the spectrum of A. The spectrum of L on F consists of five bands: 
[-0.830,-0.606], [-0.518,-0.297], [-0.219,0.219], [0,0.485], [0.611,1.000], with the 
third and fourth bands overlapping. The second branch is similar in appearance to 
the one of the operator A (see Figure H]). It shares the property that the band edges 
occur away from the set S. 



-0.297 



-0.518 



Figure 4. A gray scale image of the second branch for the operator L. Extrema are 
highlighted. 



Analogously to the case of the operator A, we again find that the extremal values 
on S of the second band function can only occur at the points k ^ X. It turns out 
that the maximum value on S is —1/3 ~ —0.333 at k = (0,0) and k = (0, vr), and 
the minimum value is — -\/2/3 ~ —0.471 at k = (vr, 0) and k = (7r,7r). These can be 
compared with the numerically found maximum over all of B of approximately —0.297 
at ^ (5.40,0.88) and -0.518 at A; ^ (0.26,0.88). The difference is 10.8% at the upper 
edge, and 9.5% at the lower edge. 

In summary, we have described two difference operators A and L acting on a 
periodic graph F, which have spectra with band edges occurring away from the boundary 
of the Brillouin zone. 
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2.3. An example in the presence of point symmetries 

The previously described examples dealt with a graph F that had only translation 
invariance with respect to Z^, and no point symmetries (i.e., symmetries that would fix 
a point on the graph). Thus, the Brillouin zone was not reduced. In this section, we 
provide an example where the point symmetry group is non-trivial, while the effect we 
observed in the previous sections still holds. We will also observe in some cases that 
spectral edges can occur on S, but not on X. 




Figure [5] depicts the fundamental domain of a more symmetric periodic graph. 

We now define the Floquet Laplacian A{k) : i'^{V) i'^{V) as the restriction of A 
defined as in (12.11) to the space il- In terms of the basis of the delta functions at vertices 
of V, this operator has the following matrix: 



A{k) :- 



where A is the matrix 




A{k) :- 



(2.10) 








1 




1 
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(2.11) 
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and B is 



( 1 











o\ 













































^0 











1/ 



(2.12) 



The matrix B describes the interaction between the two symmetric halves of the 
fundamental domain V . 

In a similar way, one can define Lilz) and observe that 



where T is the matrix 



T- 








(2.13) 



/2 











\ 






















75 

















2 
















V3 ) 



(2.14) 



The fundamental domain has 10 vertices, and so there are 10 bands to the 
spectrum. The lowest 5 bands are approximately [—3.840,-2.265], [-2.943,-1.834] 
[-1.865,-1.113], [-1.536,-0.333] and [-0.803,0.377]. The spectrum is symmetric 
about A = 0, and the remaining five bands are reflections of the previously mentioned 
bands about this point. 

We again focus on a single example: the third band [—1.865, —1.113]. A greyscale 
plot of the reduced Brillouin zonelJi] for the solution curve corresponding to this band is 
given in Figure El On the other hand. Figure [7] represents the band [-2.943,-1.834], 
for which maxima and minima occur both on the boundary H of the reduced Brillouin 
zone (albeit, not on X). For the graph A we have observed that only the lower edge of 
the third band, and upper edge of the eighth band are away from the boundary of the 
reduced Brillouin zone. All other band edges do occur on these lines of symmetry of the 
Brillouin zone. 

In order to check that the minimum point of the third branch really occurs away 
from the boundary of the reduced Brillouin zone, we repeat the procedure described in 
Section [2m With the aid of symbolic algebra computer packages such as Maple it is easy 
to find the characteristic polynomial of the 10 x 10 matrix f l2.10p . and to compute the 
derivative along the lines fci, = 0, vr. This leads to a set of four pairs of polynomial 
equations (in the variables A and cos(/i;i) (or cos(A;2))) to be solved simultaneously. 
Numerical root finding of this system reveals the minimum along the boundary to occur 
at the point k ^ (1.970,0) attaining a minimum value A ~ —1.830. We compare this 

II Note that for this figure, the reduced BriUouin zone [0, tt]^ is plotted. The picture for the full Brillouin 
zone is obtained by reflection of this picture in 2 directions. 
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-1.113 



-1.865 



Figure 6. A gray scale image of ttie ttiird brancli of the spectrum of A(fc). Extreme 
points are highlighted. 




-1.834 



-2.943 



Figure 7. A gray scale image of the second branch of the spectrum of A(fc). Extreme 
points are highlighted. 



with the strictly lower point in the interior at A; ~ (1-41, 1.78) which takes the value 
A ~ —1.865, a difference of approximately 2% (or 5% of the width of the band). 

Looking only at the four "corner" points: (0,0), (0,7r), (7r,0), {it, it), the minimum 
value attained on the third band is A ^ —1.568. 

For the operator L{k) the situation is qualitatively similar. The positions of 
the various maxima and minima do not move much. The minimum of the third 
band is located again away from the boundary of the Brillouin zone. At the point 
k ^ (1.28, 1.93) the value attained by the third band function is A ~ —0.5486. It turns 
out that the minimum value taken on the boundary of the reduced Brillouin zone is 
A ^ -0.5380 at A; ^ (1.92, 0), a difference of 2% (or 6% of the width of the band). 

One can also notice that we encounter here examples of spectral edges occurring 
on H, but not on X. So, all the possibilities do materialize: spectral edges occurring on 
X only, on S — X, and finally on B — S. 
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Figure 8. A graph of the third branch. 



3. Quantum graph case 

In this section we consider the spectrum of a periodic quantum graph G with the same 
topology as F. The spectrum a{G) can be related to the spectrum of the Floquet 
Laplacian L{k) investigated in the previous section. As a consequence, we will discover 
that the maxima and minima of some branches, and thus spectral edges as well, occur 
at the same quasi-momenta in both systems. Hence, spectral edges of this periodic 
quantum graph Hamiltonian occur inside the Brillouin zone. 

We construct a metric graph G by equipping all edges of F with unit lengths. To 
complete the definition of a quantum graph, we need to define a self-adjoint differential 
Hamiltonian. As such, we consider the negative second derivative on the edges e of G: 




(3.15) 



The Hilbert space where the operator acts is L'^{G) — 0gL^(e). The domain of 
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the operator consists of all functions / such that 
/ G H^{e) for each edge e, 

(3.16) 

/ is continuous at each vertex 
^^/g(f) = (KirchhofT, or Neumann, conditions). 

Here, /e(f ) denotes the outgoing derivative of / at f along the edge e. 

It is well known (e.g., |T51 HOI SSI SSI ) ^^^1 easy to check that Floquet theory applies 
to the quantum graph case. In particular, the spectrum ct{H) coincides with the union 
over the Brillouin zone B of the spectra of Floquet Hamiltonians H{k). Here H{k) is 
the operator defined similarly to H on Hf^^ functions with the same Kirchhoff vertex 
conditions as in (13.161) . and with the additional cyclic (Bloch, Floquet) condition (12. 3p : 

f{T{px)) = e''-Pf{x) (3.17) 

for any p G Z^, x E G. We will call such functions Bloch (generalized) eigenfunctions. 
Thus, describing the spectrum of an either combinatorial, or quantum periodic graph 
operator, we can work with such generalized eigenfunctions only. 

So, let be a Bloch eigenfunction of if on G with a quasimomentum k and 
eigenvalue lo"^, 

Hip = uj^ij . (3.18) 

Let us define a function (p on the combinatorial counterpart F of G as the restriction 
of ip to the vertices of G. Clearly, (p is a Bloch function with the same quasimomentum 
k (e.g., [IQ]). Due to (I3.18p . on each edge e = {u,v) the function ip can be written in 
terms of its values (p{v), (p{u) at the endpoints: 

/ ip(u) — cosuj ip(v)\ 

We[Xe) = '^[Vj COSUJXe + I : I SUlUJXe (3.19) 

Here and Xe denote the restriction of ip to the edge e and the coordinate along 
e respectively. Then one has 

^',{v) = (3.20) 
tancu 

Vertex conditions (13.161) imply that at each vertex v the equation 

(u) = COSCo'(y9(f) (3.21) 



LL"-' U 

holds. Thus, (f is a Bloch eigenfunction of the difference operator j- fi'^) T 
with eigenvalue cos a;. Notice that the spectrum of this operator coincides with the 
spectrum of its symmetrized version L investigated in the previous section. Thus, we 
have constructed a quantum graph example of a periodic operator with a spectral edge 
attained inside the Brillouin zone. 
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The reader might notice that the correspondence between the spectra of L and H 
works smoothly only outside zeros of sinu;, i.e. not on the Dirichlet spectrum of the 
edges. This is a well known phenomenon, see for instance [ID]. However, it is easy 
to observe that this does not influence our case and thus we have an example when a 
spectral edge of a periodic quantum graph occurs inside the Brillouin zone. 

4. Neumann waveguides and periodic tubular manifolds 

In this section, we will show existence of periodic elliptic second order operators on 
manifolds with a free co-compact action of Z^, some of whose spectral edges are attained 
inside the Brillouin zone. The simplest example is of the Laplace operator with Neumann 
boundary conditions in a periodic planar waveguide. 

In order to construct the guide, let us assume our graph G (see Figure [1]) to be 
embedded into the plane in such a way that: 

1. Each edge is a smooth simple curve of length 1. 

2. Edges intersect only at the vertices. 

3. Edges intersect transversally (i.e., there are no tangent edges). 

4. The embedded graph is Z^-periodic. 
Such an embedding is clearly possible. 

Let us take a small e > and consider a "fattened graph" domain that consists 
of tubular neighborhoods of the edges (domain U in Figure [9] below) and neighborhoods 
of vertices (domain V in Figure [9]). 



We will assume the following conditions on the domain Q,^: 

1. The boundary is sufficiently smooth (e.g., C^). 

2. The domains U have constant width e in directions normal to the edges. 

3. The vertex neighborhoods V satisfy the following property: there exist balls 6^ and 

of radii re and Re correspondingly, centered at each vertex and such that <Z V G B^. 
Besides, V must be star-shaped with respect to all points of b^- 

4. The domain is Z^-periodic. 




Figure 9. A "fattened graph" domain. 
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It is easy to see that one can construct an e-dependent family of domains satisfying all 
these properties. 

Consider now the (positive) Laplace operator —A jv^e in fi^ with Neumann boundary 
conditions on dfl,.. 

It is proven in [HlllS], [5l]-[57j that for any value of the quasimomentum k and any 
finite interval I of the spectral axis, the part in / of the spectrum of the Floquet operator 
—Afyf^^{k) converges to the corresponding part of the spectrum of the quantum graph 
Hamiltonian ( I3.15p - p.l6p on G. Moreover, this convergence is uniform with respect to 
k. This, in particular, implies immediately the following 

Theorem 4.1. For sufficiently small values of e, there is an isolated band of the 
spectrum of the waveguide operator — Ajv,e, whose end points are attained strictly inside 
the first Brillouin zone. 

Proof. Indeed, this property holds for the quantum graph Hamiltonian H on G, and 
thus the convergence result shows that it survives in fi^ for small values of e. □ 

This construction does not necessarily require the graph to be planar. For instance, 
one would not be able to have a planar embedding with required properties for the more 
symmetric graph A considered in Section 12.31 However, one can embed A into in 
such a way that it is Z^-periodic, with all other properties as required before. Then a 
3D waveguide domain can be constructed around A in a similar manner to the one 
above, such that the statement of Theorem 14.11 still holds. 

Another type of examples can be constructed as a "tight sleeve Riemannian 
manifold" around graphs G or A. The notion of such a manifold can be easily 
understood from the Figure [TO] (see precise definitions in [17] ). 



Such a manifold can be constructed preserving a co-compact free action of Z^. Then 
the results of [I7j concerning convergence of spectra of the Laplace-Beltrami operator 
on Me to those on the graph, show that the following result holds: 

Theorem 4.2. In any dimension d ^ 2, there exists an example of a closed d- 
dimensional manifold M with a co- compact, free, isometric action of 1? , such that 
there is an isolated band of the spectrum of the Laplace-Beltrami operator —Am, whose 
end points are attained strictly inside the first Brillouin zone. 





Figure 10. A "sleeve" manifold. 
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The proof coincides with the one of the previous theorem. 
5. Schrodinger and Maxwell operators 

The previous discussion does not leave any doubt that examples can be found for 
essentially any type of periodic equations of mathematical physics. However, it is 
desirable to have such explicitly described examples for the cases of periodic Schrodinger 
and Maxwell equations, interest in which stems from the solid state and photonic crystal 
theories (e.g., [1], [I9]-[2T], [301 [31], [35]-[37], [53]). 

Although we do not currently have rigorous arguments to show the existence of 
such examples, we can expect that they may be obtained as follows. Consider a planar 
embedding of the graph G, as the one considered in the previous section, with the 
additional requirement that at each vertex the tangent lines to the converging edges 
form equal angles. This can obviously be achieved. Consider then a "fattened graph" 
domain described before and the Schrodinger operator S := —A + V{x) in with 
the Z^-periodic potential V{x) that is equal to zero in and equal to a large constant 
C outside. 

Conjecture 5.1. Under appropriate asymptotics e —>■ 0,C —>■ oo, the spectrum of 
the operator S will display an isolated spectral band with its edges attained inside the 
Brillouin zone. 

What is lacking here, in spite of significant attention paid to such asymptotics 
(e.g., [ni[ia[l5lEl[T8l[22l[23l[25l[36l[371E8l M-M. 
[63]). is a spectral convergence result analogous to the one for the Neumann Laplace 
operator. Moreover, it is known that such convergence (even after appropriate spectral 
re-scaling), does not hold, due to the appearance of low energy (below the energy of the 
first transversal eigenfunction) bound states attached to vertices [121 [13 EH] . However, 
for creating an example that we are looking for, the full spectral convergence is not 
truly needed. What is required, is some kind of convergence above the energy of the 
first transversal eigenfunction, which must hold (see some results in this direction in 
[m E]). When the angles formed by the edges are equal, we expect vertex conditions 
of Kirchhoff type to arise. 

Concerning the periodic Maxwell operators V x e~^(x) Vx , where e{x) is the electric 
permeability, we expect that the simplest to come by will be an example of a 2D periodic 
medium (i.e., a medium which is periodic in two directions and homogeneous in the 
third one). As it is well known [29l [30l E], [37], in this case the Maxwell operator splits 
(according to two polarizations) into the direct sum of two scalar operators — V-£~^(a;)V 
and —e~^{x)A. We expect that similar high-contrast narrow media as described above 
should provide necessary examples (see also considerations of such high contrast limits 
in [3], [l9]-[2T],[2a[2Hl[37l[lIl[l2l[60]). 
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6. Why do spectral edges often appear at the symmetric points of the 
Brillouin zone? 

In this section, we will show why, in spite of the examples of this paper, the spectral 
edges are often attained at the highest symmetry points, and hence at the boundary of 
the reduced Brillouin zone. Namely, let Hq be a periodic self-adjoint elliptic Hamiltonian 
with real coefficients (i.e., the corresponding non-stationary Schrodinger equation has 
time-reversal symmetry). Suppose that the spectral edges of Hq are attained at 
symmetry points of the Brillouin zone B = [— 7r,7r]"' only (see the details below). Then 
we will show that for a "generic" Hq, this feature of the spectrum cannot be destroyed by 
small perturbations (with the same symmetry) of the operator. This robustness might 
be the reason why one very rarely observes spectral edges appearing inside the reduced 
Brillouin zone. 

Let us now introduce some notions. We will assume, for simplicity of presentation, 
that the Hamiltonian is the Schrodinger operator in M": 



where V{x) is a real- valued bounded potential such that V{x + p) = V{x) for all 
integer vectors p G Z". For any quasimomentum k E B, we will denote by H{k) 
the Bloch Hamiltonian defined on Z"-periodic functions (i.e., on functions on the torus 

= M"/Z") as 



It depends polynomially, and thus analytically, on k. We denote by Xj{k),j = 1,2,..., 
the eigenvalues of H{k) counted with their multiplicity in non- decreasing order. The 
band functions Xj{ ■ ) are continuous functions of k E B. Since the potential V is 
real- valued, the eigenvalues are also even in k, i.e. Xj{—k) = Xj{k). This follows from 
the fact that complex conjugate to an eigenfunction is an eigenfunction (presence of a 
magnetic potential would destroy this symmetry). This symmetry will be crucial for 
what follows. 
The ranges 



are closed finite intervals of the spectral axis {spectral bands), whose union is the 
spectrum (j(ifo)- Global maxima and minima of the band functions Xj{ ■ ) are the 
endpoints (edges) of spectral bands. It is also known (e.g., [321 [35l [53]) that X/s are 
analytic in k away from the eigenvalue crossing points. In case of the crossing, we point 
out the following elementary, but useful result: 

Lemma 6.1. Let us fix an open interval A = (a, b) C M. Suppose that for a neighborhood 
U G B the band functions Xg, j ^ s ^ j + m, satisfy 



— ^ 



A + V{x), 



(6.22) 




A, = {X,{k)\keB} 



Xs{k) e A, k eU, 



Spectral edges of periodic operators 



18 



and that the remaining band functions take values in U that lie outside a neighborhood 
of the closed interval A. Then the functions 

j+m j+'fn 

W \s{k) and ^ Xs{k) (6.23) 

s=j s=j 

are analytic with respect to k E U . 

Proof. Let us assume that m = 1 (the case of arbitrary m works out exactly same way) , 
i.e. we have two eigenvalue branches A_(/c) := \j{k) and \+{k) := Aj+i(fc). Consider a 
positively oriented circle F C C centered at (a + 6)/2 with radius {b — a)/2 + e with a 
small e > 0. The two-dimensional projection 



(6.24) 



r 



depends analytically on k E U. Let M{k) be the range of P{k). It forms an analytic 
two-dimensional vector-bundle [35l |62]. Let {61,62} be a basis in M(/co) with some 
fco G U . For k close to k^ we can define a basis of M{k) analytically depending on k as 
follows: 

f,{k):=P{k)e,. (6.25) 
In this basis, the operator function 

P{k)H{k)P{k)\Mik) = H{k)P{k)\M(k) 

can be written as an analytic 2x2 matrix-function A{k) with eigenvalue branches X-{k) 
and X+{k). Therefore, the functions 

detA(fc), ti A{k) 

are analytic in /c in a neighborhood of /cq. □ 

The only fixed points k E B for the symmetries k —k + p,p E 27rZ", are the 
ones from the set 

X = {k= {k,, ...,kn)eB\k^e {0, tt}, J = 1, ■ ■ ■ , n}. (6.26) 

In view of the symmetry Xj{k) = Xj{—k) we also have Xj{kQ + k) = Xj{ko — k) for 
any fco G X. We have already shown in this text that the global extrema of the band 
functions can occur outside the set X. The experimental observation, however, is that 
for most periodic operators of practical importance and for practical values of their 
parameters (e.g., potentials, electric permittivity, etc.), the band endpoints do occur 
on X. One can easily observe this by looking at dispersion curve calculations in solid 
state physics or photonic crystals literature (e.g., |51i ?]). Our main question now 
is: Why do the spectral edges occur so often on X? 

Our considerations will be local on the spectrum. Thus, let us fix a finite interval 
A = (a, b) of the spectral axis. Note that the number of spectral bands Aj overlapping 
with A is finite. We first introduce the following notion: 



Spectral edges of periodic operators 



19 



Definition 6.2. We call a periodic Hamiltonian H simple on a finite interval A, if 
the global extrema of the hand functions Xj{k) which occur inside A, are attained at the 
points of the set X only. 

The simplicity property defined above will be discussed for "generic" periodic 
operators: 

Definition 6.3. We call a periodic Hamiltonian H generic on a finite interval A, if 
for every band edge Aq occuring inside A, the band functions Xj assume the value Aq at 
finitely many points of the Brillouin zone B, and in a neighborhood U of each such point 
ko one of the following two conditions is satisfied: 

(i) There is a unique band function A( ■ ) for k E U such that A(A;o) = Aq; moreover, 
kf) is a non- degenerate extremum of X{k). 

(a) For k E U there are only two band functions A+( ■ ), A_(- ) such that A_(/co) = 
A+(/co) = Aq. Moreover, A_(fc) < X+{k) for all k E U — {k^}, and ko is a non- 
degenerate maximum of the product D(k) = (A+(A;) — Ao)(A_(A;) — Aq). 

Above by a non-degenerate extremum we understand an extremum with a non-degenerate 
Hessian. 

Recall that in view of Lemma 16.11 the determinant of A{k) is analytic on U. 
Definition 16.31 means that the band functions of a generic Hamiltonian behave near 
band edges as eigenvalues of a "generic" 2x2 self-adjoint analytic matrix function. We 
refer to cases (1) and (2) in the above definition as the single edge case and the case of 
two touching bands respectively. Note also that in Definition 16.31 the band edge Aq is 
not assumed to be (albeit could be) an endpoint of the spectrum. 

The following conjecture is believed to hold (see a variety of similar genericity 
conjectures in, e.g., [21 El [371 HE] ) . 

Conjecture 6.4. Generic periodic Hamiltonians form a set of second Baire category in 
a suitable class of periodic operators. 

The closest to the proof of this conjecture is the result of [33], where it was shown 
that generically a band edge, which is an endpoint of the spectrum, is attained by a 
single band function. 

Our aim is to show that a generic simple Hamiltonian Hq remains generic and 
simple under small perturbations. More precisely, we introduce the family of operators 

Hg = Ho + gV{x), Ho = -A + Vo, 

where Vq and V are bounded real-valued Z"-periodic functions, and (7 G M is a parameter. 
We denote by Xj{k,g) the band functions of Hg. If g = 0, we drop g and write simply 
Xj{k). Since Hg is analytic in g, the band functions Xj{k,g) are analytic in {k,g) away 
from the crossing points, and the quantities defined in fl6.23p are analytic in {k, g) under 
the conditions of Lemma 16.11 

We can now formulate a result that gives a partial answer to the question posed in 
this Section. 



Spectral edges of periodic operators 



20 



Theorem 6.5. Let A C M 6e a finite closed interval, and let the operator Hq (see ( 16'. 2^) ) 
be simple and generic in a neighborhood of A. Then, for sufficiently small values of g, 
the operator Hg is also simple and generic in a neighborhood of A. 

Proof. Let A' be a finite closed interval containing A in its interior and such that 
operator Hq is simple and generic in an open neighborhood A" of A'. The continuity 
of \j{k,g) in g guarantees that for small g, the spectral band edges of the perturbed 
operator occurring on A', are either perturbations of the band edges of Hq that are 
inside A", or are produced by opening a gap between two touching spectral bands of 
Ho. 

Let Ao G A" be a single band edge of Hq, or the point where two bands touch. 
Assume without loss of generality that Aq = 0. Since the unperturbed operator Hq 
is simple, again, by continuity of Xj{k,g) in g, for sufficiently small values of g, the 
perturbed eigenvalues cannot reach their global maxima outside a neighborhood of the 
set X. Thus, it suffices to consider the neighborhood of each point & X individually. 
We assume without loss of generality that fco = 0. 

Further proof requires different arguments for the two cases featuring in Definition 

El 

6.1. The single edge case 

Let X{k) be the unique band function of the operator Hq, which attains at /cq = its 
non-degenerate extremum, which for definiteness will be assumed to be a maximum. 
Recall that A( ■ ) is analytic in k and X{k) = X{—k), so that 

X{k) = X2{k) + X,{k), 

where A2 is a negative definite quadratic form and Ag is an analytic function such that 
Ae(/c) = Odfcl"^). For sufficiently small g and k, the eigenvalue X{k,g) will remain 
separated from the rest of the spectrum of Hg{k). Thus, to complete the proof, we 
need to show that X{ ■ ,g) attains its maximal value at A; = and this maximum is 
non-degenerate. Due to analyticity, 

X{k,g) = X2{k) + X,{k)+gX{k,g), 

where X{k,g) is a real-valued real-analytic function of {k,g), and X{k,g) = X{—k,g). 
The latter property implies that 

VkXik,g) = Oi\k\), 

uniformly in g. Making an appropriate linear change of variables, we can always assume 
that X2{k) = —\k\^/2. Then, taking the gradient with respect to k, we obtain 

VfcA(A;, g) = -k + VkXeik) + gVkX{k, g). 

Consequently 

\VkXe{k) + gVkX{k, g)\ ^ C(|A;|3 + g\k\), 
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with some positive constant C, and hence, for \g\ < (4C) ^, < (2v^) \ 7^ 0, we 
get, 

\V,X{k,g)\>^-^^0. 

This proves that the only stationary point of X{ ■ , g) is k = 0. Moreover, since ■ ) is 
negative definite, the function A( • , g) has a non-degenerate Hessian if g is sufficiently 
small. Thus the band function X{k,g) attains its extremum on X and satisfies the 
requirements of Definition 16.3( 1). 

6.2. The case of two touching bands 

Assume, as above, that Aq = 0, k^ = 0, and that A_(fc) and A_|_(/c) are the band functions 
as given in Definition 16. 3[ Denote by X±{k,g) the perturbed band functions. According 
to Lemma [6.11 the functions 



d{k,g) = X_{k,g)X+{k,g), t{k,g) = ^(A_(A;,^) + X+{k,g)) 

are analytic in a neighborhood of {k, g) = (0, 0). Remembering the central symmetry of 
the eigenvalues and the genericity assumption for Hq, we can write 

d{k, g) = d2{k) + de{k) + gd{k, g), 

^ (0.27) 
t{k,g) = t2{k)+teik)+gt{k,g). 

Here all functions are analytic near {k, g) = (0, 0) and even in k. The functions ^2 
and d2 are quadratic forms, the terms de,te are 0(1^1"^), and, by virtue of genericity, 
d2 is negative definite. Thus, as in the first part of the proof, we may assume that 
d2{k) = -\k\^/2. Note also, that d{0,g) = 0, since the eigenvalues Xj-{0,g) are of order 
0{g), and hence d{0,g) = 0{g'^). Introduce the quantity 

m = t'^-d = hx+-X_Y ^0. 
Using (16.271) we get 

"^(^, 9) = g'^moig) - d2ik) + rrieik, g), m^ik, g) = 0{\k\'^){k'^ + g), 

where the functions niQ^g) = i{0,g)'^ — dgd{0,g) and me{k,g) are analytic in k,g, and 
me is even in k. Since m{0,g) = g'^mQ{g), we also have mo{g) ^ 0. 

Let us list some simple estimates that these functions and their gradients with 
respect to k satisfy in a neighborhood of (0, 0). Below we denote by C, Ci some positive 
constants whose precise value is not important: 

\tik,g)\^Cie + g), \mik,g)\^Ci\k\' + g'), 
\Vkt{k,g)l\Vkm{k,g)\ ^ C\kl 



(6.28) 



\Vkm{k,g)\ ^ hk\, 

^ ^ (6.29) 

■m{k,g) ^ g'^mo{g) + 
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The eigenvalues \±{k,g) solve the characteristic equation 
X"^ -2t{k,g)X + d{k,g) = 0, 

and thus 

A±(fc, g) = t{k, g) ± ^m{k,g). (6.30) 

By (16.291) the eigenvalues X±{k, g) can coincide only at A; = 0, so that they are analytic in 
k for A; 7^ 0. Let us prove that X±{ ■ , g) have no stationary points if 7^ 0. Differentiate: 

V,X^{Kg) = V,t{k,g)±^^^^P^ 
Now the estimates (I6.28P and (16.291) imply: 



- |Vfct| ^ ^ - C\k\ ^ Ci\k\ 



V\W+7 

for small g and A; 7^ 0. This proves that A-|-( ■ ,g) attain their extrema only at A; = 0. 

It remains to show that the eigenvalues satisfy the requirements either of Part (1) 
or Part (2) of Definition 16. 3[ If mo{g) > 0, then by (I6.29P and (I6.30p . the eigenvalues 
X+{k,g) and X-{k,g) are decoupled for all k and g, and their extrema are clearly non- 
degenerate. If mo{g) = 0, then by (I6.30p A+(0,5') = X^{0,g), and then by (I6.27P the 
determinant d{k,g) has a non-degenerate Hessian for small k,g. 

The proof of the 16.51 is complete. □ 



7. Final remarks 



• Suppose that for a particular periodic operator the spectral edges do occur at the 
point of X = {k\ kj = ±7r or 0} only. This means then that one can find the correct 
spectral edges (and thus the spectrum as a set), computing only spectra of problems 
that are periodic or anti-periodic with respect to each variable (say, periodic with 
respect to Xi and X3 and anti-periodic with respect to X2). This resembles then the 
ID situation [131 ES] , when the edges of the spectrum are attained at the spectra 
of the periodic and anti-periodic problems. 

• In the last Section we have restricted ourselves to the case of Schrodinger operators 
with electric potentials only. However, the proof in fact does not use the structure 
of the operator and could be extended to arbitrary analytically fibered operator in 
the sense of [26], as long as the central symmetry k ^ —k holds. 

We also assumed parametric perturbation (i.e., perturbation by gV with a small 
scalar parameter g). However, without any change in the proof, one can consider 
V a.s a. functional parameter and prove the same statements for small values of this 
parameter. 

• Observations of stability under small perturbations of critical points of a function 
in "general position" in presence of symmetries, analogous to the ones in the last 
section, have been made before in different circumstances, e.g. in |27j . 
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